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Abstract 

Extending the classic works of Berg and Purcell on the biophysics of bacte- 
rial chemotaxis, we find the optimal chemotactic strategy for the peritrichous 
bacterium E. Coli in the high and low signal to noise ratio limits. The opti- 
mal strategy depends on properties of the environment and properties of the 
individual bacterium and is therefore highly adaptive. We review experiments 
relevant to testing both the form of the proposed strategy and its adaptability 
and propose extensions of them which could test the limits of the adaptability 
in this simplest sensory processing system. 
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A. Introduction 



If placed in an inhomogeneous solution of a chemoattractant such as a-methyl-D,L- 
aspartate, E. Coli collect visibly in the regions of high concentration of the attractant 
There is a parallel phenomenon that occurs for chemorepellants where the bacteria collect 
in regions of low repellant concentration. The two phenomena are referred to as chemotaxis 
and have been known in a variety of bacteria since the 1880's ||. They constitute perhaps 
the simplest known example of sensory processing dependent behavior in a living organism. 
Actual sensory processing must be involved since, for E. Coli, the chemoattractant need 
not be a substance that the bacterium can metabolize in any way and it is known that the 
response to the chemoattractant relies on specific chemoreceptors on the outer membrane 
of the bacterium which bind chemoattractant and signal their occupancy to the interior 
of the cell through a phosphorylation cascade |L[||. The manner in which the bacteria 
use this information to reach the regions of high chemoattractant concentration was first 
illuminated by Berg and Brown || who built a special microscope to track the motions of 
the individual bacteria. What they saw were stretches of motion at approximately 10/im/sec 
with a slowly varying direction of orientation (termed "runs" ) separated by periods when the 
bacterium came to a stop and changed orientation very rapidly (referred to as "tumbles") 
before continuing on in another run. We now know that both of these characteristic motions 
of E. Coli, and ultimately chemotaxis, are due to the rotation of the flagella [g-§|. E. 
Coli typically have several flagella spread over their surfaces (this is what is meant by the 
designation peritrichous) and runs result from counterclockwise rotation of all the flagella. 
For that direction of rotation, the flagella come together to form a "bundle" and cooperate 
to propel the bacterium. Because the flagella are helical, the opposite sense of rotation has 
quite a different effect: the flagellar bundle comes apart, the bacterium is not propelled and 
its orientation varies rapidly, but apparently randomly ||, resulting in a tumble. Chemotaxis 
results from the coordination of the tumbling times with the time dependence of the receptor 
occupancies so that the bacteria change direction less often when they are headed in the 
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direction of increasing chemoattractant |T0|,|9[|. The problem we will discuss in this work is 
the optimal strategy for coordinating these tumbles with the input from the chemoreceptors. 

Clearly, for an optimal strategy to exist at all a problem must be very highly constrained. 
The principle constraints which make this problem solvable are taken largely from the exper- 
imental literature on chemotaxis and motility of E. Coli. First, as was pointed out by Berg 



and Purcell |TT[ in this context and as we will briefly discuss, due to rotational Brownian 
motion, E. Coli can not maintain an orientation for an extended period of time. Second, 
E. Coli make no controlled changes of direction. It is obvious, given that they can not 
maintain orientation that, for sensory processing reasons alone, they are incapable of turn- 
ing in a specific direction |lT|,|5|, however a change of direction of a controlled magnitude 
is in principle possible. In practice, there is some evidence that the length of tumbles is 



affected by sensory input under some circumstances |12| ; however, this is not believed to be 
important for chemotaxis under realistic conditions ||. We therefore make the assumption 
that E. Coli change direction by entering into "tumbles" which have have no characteristics 
which depend on sensory input. It seems likely that, in view of the limited use E. Coli could 
make use of steering capability given its orientation problems, this simple method of direc- 
tion change was evolutionarily preferred because the cost associated with this capability are 
lower than those that a more developed steering capability would impose. In any case, we 
assume the tumbles are all identical and effectively randomize the orientation of the bacteria 
over the course of a time r tumble ~ .15 sec [[L3|]. Finally, we assume that during runs the 
bacterium swims with fixed velocity, v , independent of sensory inputs. This is known to be 
approximately true experimentally ||. 

Taken together, these three constraints are sufficient to allow us to determine the optimal 
chemotactic strategy in the high signal to noise ratio limit for the following definition of 
optimality: we call a chemotactic strategy optimal if it maximizes the expectation value of 
(v ■ Vc) for a static, uniform gradient. Here v is the swimming direction of the bacterium 
and the average is over the time history of the bacterium's path. There are alternative 
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definitions of optimality such as that of maximizing (c) or (here the r subscript denotes 

an average only over the runs ||14||), but we expect them to result in very similar optimal 
strategies provided they do not incorporate directly (1) game-theoretic competition between 
bacteria or (2) a complex structure of maxima and minima in the concentration so that 
gradient descent approaches such as we are proposing become trapped in local minima. The 
neglect of the latter possibility seems very reasonable for realistic environments, however, 
considerations like the former may well have played an important role in the evolution of 
E. Coli. For example, the known pattern forming behavior of E. Coli [T5|] demonstrates 
that the problem of chemotaxis has special properties in presence of depletable nutrient 
and large numbers of other E. Coli. In fact, our choice of optimality is partly motivated 
by the consideration of competition. Maximizing (v ■ Vc) for static, uniform gradients 
chooses the strategy that leads the bacterium to attractant most rapidly, which is probably 
evolutionarily preferable to one that leads the bacterium there more slowly, but then results 
in the bacterium staying slightly closer to the region of maximal concentration. 

Before discussing the strategy we obtain, let us first make a general remark on the nature 
of any optimal strategy: with the definition of optimality we have chosen, the strategy must 
consist of a deterministic algorithm for deciding when to tumble based on the history of 
chemoreceptor occupancy. For any given history, the bacterium's expected, future, average 
v ■ Vc is either raised or lowered by initiating a tumble at that particular moment and if 
it is raised the bacterium should tumble and if it is lowered it should not. The stochastic 
nature of observed runs and tumbles should result (for an optimal strategy) entirely from 
the stochastic nature of the inputs (chemoreceptor bindings) , not from any deliberate intro- 
duction of "noise" in the decision making on the part of the bacterium. In practice, a totally 
deterministic strategy at fixed input requires an amplifier of arbitrary fidelity and gain to 
allow the inputs to drive the decision making apparatus and is not realizable; however, the 
phosphorylation cascade seems capable, in practice, of providing sufficient gain and fidelity 
that the tumbling would be effectively deterministic. A deterministic strategy is, however, 
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in some conflict with the two state model proposed in [0 on experimental grounds, where 
the past history of the receptor occupancies is taken to modulate the rates with which the 
flagellar motors change their direction of rotation. In that strategy, it is not the states but 
the the rates for transitions between the states that are set for each flagellar motor and they 
are all set separately, so that the only correlations between different flagella are rate-rate 
correlations. In that case, the strategy would actually be stochastic even for fixed inputs to 
the chemoreceptors. However, in the limit that the rate modulations for going from running 
to tumbling are large (the rate is either ~ or ~ oo) the two state model has a strategy of 
the form we propose. 

As mentioned above, in practice, a purely "deterministic" strategy for switching between 
two states is impossible and, to approach it, large variations in the rate of the above type 
could be used. The rate model is thus one possibility for a realistic process which approaches 
the optimal strategy and may represent E. Coifs best effort to implement the optimal 
strategy. However, it should be noted that the evidence for the modulated rate model is 
not entirely conclusive. In particular, as we will see in our discussion of run and tumble 
statistics for the optimal strategies, the exponential tails present in the durations of run and 
tumble times are not uniquely explained by the rate model, as claimed in |TJ[ . Further, there 



is a pronounced advantage that a deterministic strategy has when the problem of flagellar 



coordination is considered 16 



Free swimming E. Coli in the absence of gradients spend roughly a tenth of their time 
tumbling || while tethered bacteria |20| with a single flagellum rotate it clockwise nearly 
half of the time. Clearly, if the typical bacterium has five flagella, there must be significant 
conflict between the flagella in the absence of any signal. The formation of a coherent flagellar 
bundle 80% of the time requires some sort of coordination. Notice that no coordination can 
result from the signals to the flagella if the signals are independent, as the model of (12 



predicts they should be at small signal to noise ratio, where the above numbers apply. At 
least one mechanism for the required coordination has been proposed [[11]] and it seems 
clear that some explanation is required, unless the data from tethered bacteria are taken to 



be unrepresentative. If the flagella are coordinated then much of the self-induced "noise" 
implied by the stochastic nature of their individual biases will be eliminated by the pooling 
of their inputs: an effectively deterministic strategy will result. In fact, any "deterministic" 
strategy would have to result either from some sort of cooperativity effect (to reduce the 
noise inherent in the stochastic nature of the binding and unbinding of individual internal 
signaling molecules to receptors) or from sufficient amplification of the input signal to drive 
very large variations in the concentration of internal signaling molecule. Cooperativity is 
generally a more efficient and robust solution. 

It could be supplied by a either coordination of the different flagella and pooling of their 
signals, cooperative signaling to the individual flagellum or both. In fact, the protein FliM 



which is believed to be responsible for translating the internal tumble signal |17]], coming 



from the phosphorylated form of the protein CheY |18[], is believed to be present in many 



(about 100) copies for each flagellar motor [19|. It is not known how many of these copies 



actually bind CheY or how this binding is transduced into a tumble signal, but there is 
clearly great potential for cooperativity. For example, if the tumbling/running switch were 
whether more or less than 50 molecules of phosphorylated CheY were bound, then the 
behavior would be nearly deterministic for even moderate variations in the concentration of 
phosphorylated CheY. 

In view of this, it is not clear that there is much difference in practice between a de- 
terministic signal and stochastic signals to the individual motors or, perhaps, to the many 
individual CheY binding cites at each motor, since the results of pooling the stochastic 
signals might be nearly deterministic. In so far as there is a difference, the deterministic 
strategy will lead to better performance, but may be an expensive capacity for E. Coli to 
maintain. 

We note that the "response regulator" model proposed in pl| , in which the tumbles 
are induced by threshold crossing some functional of chemoreceptor binding histories, is 
of the appropriately deterministic type and our optimal strategy will be a realization of 
this strategy where we specify, at high signal to noise ratio, the correct threshold (with the 
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functional being somewhat arbitrary). At low signal to noise ratio, both the correct threshold 
and functional of the receptor histories will be determined, subject to certain assumptions 
about the receptor correlations. We will see that the statistics of runs and tumbles resulting 
in both cases are not inconsistent with the statistics observed experimentally, contrary to 
the claim of |T2| that threshold crossing strategies, as opposed to rate modulation strategies, 
disagree with the data. 



B. High Signal to Noise Ratio 

Let us now begin with the "high" signal to noise ratio case, defined by the bacterium 
being able to measure the projection of the gradient onto its swimming direction in a time 
much shorter than any other relevant time scale, in particular the time scale set by rotational 
diffusion. For an object undergoing rotational diffusion 

(h(t) • n(0)> = exp{-2D rot t) (1) 

where D TOt is the rotational diffusion constant. For E. Coli it is know empirically to be 
about .15 radian 2 /sec, implying a time scale for rotational Brownian motion of about 3 sec. 
This is in rough agreement with what one expects for the case of a sphere of radius r = 1/zm 
in room-temperature water. There 

D rot ~ kT/ u TOt (2) 

u TOt ~ 87rr/r 3 (3) 

~ 2.5 10" 13 cm 2 /sec 

-D rot ~ -16 radian 2 /sec (4) 

where T is the temperature, z/ rot is the rotational drag on the sphere, and rj is the viscosity 
of water. Note that, since we are at low Reynolds number, the rotational diffusion which is 
disorienting the bacterium is Markovian. Its present orientation embodies all of its knowledge 
about the future and there is no need to keep track of an angular momentum. Because all 
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of the information about the future is contained in the present orientation, or equivalently 
the bacterium's present knowledge of v ■ Vc obtained from the chemoreceptor histories, the 
optimal strategy in the high signal to noise ratio limit requires only the specification of a 
single number: a threshold, T, which is the maximum value of v ■ Vc which the bacterium 
will tolerate before tumbling. Any processing strategy for the chemoreceptor histories that 
allows the determination of v ■ Vc rapidly compared to other time scales is acceptable. This 
makes this limit particularly simple to discuss. 

To solve the high signal to noise ratio limit we need to make the assumption that each 
tumble is perfectly disorienting and totally randomizes the bacterium orientation. This is 
approximately true and allows a full solution of this case. This is interesting even though 
the bacterium may rarely, if ever, encounter a high signal to noise ratio environment because 
the optimal strategy can be solved for completely and it gives us an idea what the low signal 
to noise ratio strategy is evolving towards as the signal to noise ratio is increased. In this 
limit, there is still an interesting strategy because the bacterium can not steer, only reorient 
itself through stereotyped tumbles, each of which requires a finite amount of time. This 
finite time cost will be very important in determining the optimal strategy and will result 
in the bacterium displaying a surprisingly large amount of "optimism", by which we mean 
that bacteria which know they are not swimming directly up the gradient should choose to 
continue running because of the chance that rotational diffusion will improve their prospects 
faster than tumbling would. 

Let us consider first the case where the bacterium is in a uniform gradient of chemoat- 
tractant and therefore has had time to measure the magnitude of Vc. In this case, it can 
translate it's knowledge of v ■ Vc directly into knowledge of 0, the angle between its motion 
and the concentration gradient, and the strategy consists of an angle at which the bacterium 
"decides" it is moving sufficiently in the wrong direction that it is worth taking the time 
and the risk of orienting even further in the wrong direction which a tumble will involve. 

It is clear that the dimensionless number -D ro t r tumbie fixes 6 C and one would expect that, 
for small values of -D rot r tum bi c , C would also be small, vanishing like (i} ro t ^tumble ) p - In 
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practice, -D rot r tum bi c ~ (.15sec _1 )(.15sec) ~ .02, which is indeed small, and we might expect 
that 6 C would also be small. However we will see that this is not the case as p = 1/6 and 
the prefactor between (-DrotTtumbie) 1 ^ 6 and C is important. The correct value of 6 C can be 
obtained from the following argument: the equilibrium distribution of orientation of running 
E. Coli for a specific choice of 6 C can be obtained from solving the heat equation on the 
surface of a sphere with a perfectly absorbing boundary condition at C and a source term 
that deposits new E. Coli uniformly over the allowed region at a rate that exactly cancels the 
current into the boundary at 6 C . From the distribution, one can compute (v • t^t) by first 
computing the average over the distribution of running E. Coli and then multiplying by the 
fraction of E. Coli which are running. The latter is given by the integral of the distribution 
over the allowed region of the sphere divided by the current into the boundary of the allowed 
region times the average length of time between runs. The last factor is given by the average 
length of a tumble times a factor, (1 — forbld ^" 1 arca )- 1 ) which reflects the chance for a tumble 
to result in an orientation past 6 C , in which case it will be followed immediately by another 
tumble. The relevant heat equation on the sphere reads: 



C can be fixed by the requirement that p vanish at C . The resulting, unnormalized solution 
is: 



Both C and p can be multiplied by an arbitrary, common factor since we have not yet 
imposed any normalization. The "current " into the boundary is given by 




(5) 




(6) 



C = Aot/2 



(7) 



Jloss = 27rC(l - COS0 C ) 



(8) 



while the integral of p is given by 
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2tt J " d{- cos e)p(6) = tt(cos 6 C - 1 - 4 In (««y ) ) (9) 

The expectation value of v ■ is given by 

Vc, , sin 4 % , 

f^U^Vc ^ s-T 10 

|Vc| 1 2Aotrtumbic - 2(sin 2 f + 2 In cos ^) 

In general, the maximum can be found numerically, but for the special case of -D rot r small, 

where one expects small G c , one can expand the trigonometric functions to find for x = O c /2: 

Vc . .„ . 480x 4 - 80x 6 + 6s 8 + . . . , . 

\V ■ 7=-r) ~ V VC - A 11 

X |Vc| 7 1 1 15360 Aotnumbic + 480x 4 + x 8 + .. . V ; 

anticipating the result that -Drot^tumbie oc x 6 and expanding again one finds that the derivative 
vanishes for: 



x ~ V 6D rot r tum bie (12) 



6 C ~ 2Y / 6D ro t'r t umble (13) 

(14) 

already for -D rot r tum bi c ~ .01, O c is not small (the above formula is valid only if it is small 
and would predict ~ Air already at this point). The value of C obtained from numerically 
finding the value which maximizes the return in depicted in Fig. [l|. It rises very rapidly to 
the vicinity of tt/2 where it remains for all reasonable values of -D ro tT tum bi e , until ultimately 
the behavior for very large -D rot r tum bi e is given by 



©c = 7T - \/2/-DrotTtumble (15) 

We see that, under a broad range of circumstances, the optimal strategy is essentially to 
continue if one is moving into regions of higher concentrations of attractant and to reorient 
if one is moving into regions of lower concentration of attractant. The bacterium should be 
surprisingly tolerant of moving in the directions that are far from perfectly aligned with Vc, 
even if it has perfect information as to how far off course it is. In fact, the fraction of time 
the bacterium spends tumbling, /, is not ~ 1, as one would naively expect, but rather is 
given by: 
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r _ ^tumble J loss f T ft\ 

= (1 - forbidden area ,) J 2n ^_ CQS Q) + 1 ' 

2Z) ro tTt umD i e 
cos 6 C - 1 - 4 In (cos + 2D rot T tumb i 

5 -Dxot ^tumble (18) 



(17) 



1 + 5-D ro tTt um ble 

-.1 (19) 

where we have used C ~ n/2 and DrotTtumbie ~ -02. High frequency noise added to the 
thresholded quantity will make the bacterium even more tolerant of a signal indicating that 
it is swimming the wrong direction. 

Before discussing noise, we first mention what strategy is appropriate if the bacterium is, 
for some reason, such as being in a spatially non- uniform gradient, unable to estimate the 
magnitude of the gradient and therefore the absolute angle between its present orientation 
and the gradient. First, note that, for the appropriate value of -DrotTtumbie, the critical 
angle is close to tt/2 so that, approximately, only the sign of v ■ Vc is important and the 
ignorance about |Vc| is unimportant. In so far as it is important, the optimal strategy in 
the absence of information about |Vc|, depends on the a priori probability P(|Vc|) and is 
therefore somewhat non-universal, depending strongly on the statistics of the bacterium's 
environment (it should consequently also be highly adaptive). In particular, the optimal 
strategy is now to choose a value of s CT iticai = v ■ Vc at which to tumble and the optimal 
value can be chosen by maximizing: 

(v- ]||> ~ / d(\Vc\)P(\Vc\)R (e c = arccos(^)) (20) 

where R(Q C ) is defined by the right hand side of Eq. |1(]. The quantity to be maximized 
corresponds to the average of the previously calculated return over the actual strategies 
that will result from a given choice ^critical- Clearly, as P(|Vc|) becomes sharply peaked, 
this reduces to the case where |Vc| may be taken to be known, in which case, as we have 
seen, s cr i t i C ai ~ 0; when the uncertainty in |Vc| is large, the bacterium should select an 
even smaller s cr iti ca i- The exact details depend on P(|Vc|) in the environment to which the 
bacterium is adapted. 
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Let us now discuss, in the case of a sharply peaked P(|Vc|), the introduction a small 
amount of high frequency, Gaussian noise to the bacterium's information about 0. By the 
noise being small we mean that its magnitude is much less than that of the signal v |Vc|, and 
by high frequency we mean that the characteristic time for v ■ Vc to to vary due to diffusion 
~ Q 2 /D TOt is much larger than the noise time constant, r noise , defined by 



JduN{u) 



~ \jduj uo 2 N{uo) 

Here N(u) is the power spectrum of the noise so that the mean squared magnitude of the 
noise, (r/ 2 ), is given by 

(V 2 ) = J ^N{u) (22) 

Such noise might arise inside the cell, involving the thresholding mechanism itself or the 
phosphorylation cascade, or it might come from the fluctuations in receptor occupancy if 
the cell does not completely low pass filter the receptor inputs. To compensate for the noise 
the optimal strategy for the bacterium is to set its threshold, T, on the input, v- Vc+r], where 
77 is the noise, slightly lower than in the noiseless case where without noise T = v cos C | Vc|. 
To determine T note that the rate at which v ■ + 77 crosses the threshold, T, in one 
direction is given by 

r( e) = ^ e xp(- ' r -— Q' VC " 2 ) (23) 

The ideal value of T can then be computed by converting our previous heat equation with a 
perfectly absorbing boundary into one with heat absorption at angle G given by p(Q)R(Q). 
For t>cosG|Vc| > T, R(Q) is given by r(G) and for t>cosQ|Vc| < T, by (27rr noise ) _1 . 
Treating the tumble rate as heat absorption in this way is exact if the tumbling is a mod- 
ulated Poisson process. This is true in the limit r noise — > (a sensible answer in this limit 
also requires (rj 2 ) — ► 0), but the treatment is approximate for noise with a finite correlation 
time. However, we are interested in the short correlation time limit and the approximation 
is justified. In this limit, the heat effectively can not penetrate into the region where 
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(■ycos6|Vc| — T) 2 < 2(r] 2 ) \n{D TOt r noise ) 



(24) 



equivalently 



6 > arccos 



T \ _ 2(^)| ln(27rDr noise ) 



(25) 



v\Vc\J \ v 2 \Vc\ 2 -T 2 



It is important to realize that, at this point, the tumbling rate is rapidly increasing because 
the value of the argument of the exponential, ~ | ln(27i\Dr noisc )| for small r no i S e, is large and so 
is its derivative with respect to 6. The point effectively specifies the location of an absorbing 
"wall", and, therefore, we are back to the case where our original analysis applies and the 
optimal strategy is one which chooses T so that 6 lies at the 6 C of the noiseless problem; 



this can be done by choosing the threshold to be at vcos6 c |Vc| - J 2 ^ 2 ^^^ n ° isc)l . The 



threshold has been moved further from v\Vc\ because the noise results in tumbles being 
triggered prematurely. Thus, the bacterium should be even more tolerant of an input signal 
which indicates that it is swimming in the wrong direction than we observed that it should 
be in the noiseless case. 

At this point it is worth making a few remarks about the statistics of the runs and tumbles 
in the high signal to noise ratio case. For simplicity, let us consider the noiseless case with 
known |Vc|. The bacteria emerge from tumbles with random orientations and for some 
tumbles the orientation immediately after a tumble is past the critical allowed orientation 
and they immediately tumble again. This will lead to a renormalization of the effective 
tumble length in this model from r tum bi c to r tum bi e -i — ^—et, however, the tumble durations 

_L ~ ~ COS Oc 

will continue to be purely exponentially distributed, if they were initially exponentially 
distributed, as they are found to be experimentally. Meanwhile, if we average over all the 
bacteria which do not immediately tumble to find the statistics of run duration, we find 
these are also roughly exponentially distributed. To see this, recall that the heat equation 
we solved to get p(0) is of the form: 





(26) 



where L is a linear operator and we impose the boundary condition p(0 > 6 C ) = 0. This 
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could have been solved by solving for the orthonormal eigenfunctions, /a(@), of L which 
vanish at 9 C , then re-expressing the linear source term, C, as a sum of these: 

C = Y,cxh(e) (27) 

A 

c x = / eC rf(cose)/ A (0) (28) 
Jo 

Then we have: 



p(0) = 5> A / A (0) (29) 

A 

ax = -y (30) 

The distribution of intertumble intervals is then proportional to J2x c x&~ xt - The eigenfunc- 
tions for a general C are hypergeometric functions 2 i 7 i(— f + 1; 1; C) where ( = cos G) 
p3| and solving the boundary condition is in general impossible. However, for the special 
case O c = 7r/2, the eigenfunctions are those Legendre polynomials, P ra (cos0) which vanish 
at 7r/2. The eigenvalues are Dn(n + 1), where any odd n are allowed. The smallest eigen- 
value is 2D, and the tail is therefore of the form e~ 2Dt . The next highest eigenvalue is 12.D 
and the coefficient in front of it at t = 0, C\2Di is ^ of the contribution from the n = 1 term, 
so the approach to exponential decay of the form e~ 2Dt is very rapid. This is consistent with 
the fact that the mean duration of a run is, in general, given by cos ® c 1 41n ^° s 2 which f or 

' to ' to J D rot (l-cos C ) ' 

B c = 7r/2 gives 21n ^ -1 ~ .38-D -1 . This is only slightly smaller than \D _1 because of the 
inclusion of the higher n terms. Roughly half of the shift is accounted for already by the 
inclusion of the n = 3 term. 

Although we do not believe that high signal to noise ratio limit applies to the experiments 
of |5,10,12], we should point out that the exponential tails observed for the run distribution 



is consistent with the statistics of runs and tumbles observed there. The fact that our 
procedure for ending runs is essentially of the threshold crossing type does not imply that 
the distributions of run durations has a power law tail. This contradicts the claims of |12]|, in 



which power law tails were found for the threshold crossings of a random walk and argued to 
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be typical for threshold crossing processes. A random walk is a very special case since there 
is no finite correlation time for the displacement. Threshold crossing a variable with a finite 
correlation time results in a tail for intervals between crossings that is generally exponential. 
The decay constant is roughly given by r~J relations ln(prob. no crossings in r CO rreiations), and 
only for divergent or vanishing r corre i ations should non-exponential tails result. Thus, there is 
not, in general, any a priori conflict between response regulator models and the statistics of 
runs and tumbles. 



How should the bacterium make use of the information available to it in small signal 
to noise ratio limit? As in the high signal to noise ratio case, the optimal strategy is a 
deterministic algorithm for generating tumbles based on the chemoreceptor binding histories. 
In this case, however, the bacterium will not tumble at fixed angle because it can not 
determine its orientation with respect to the gradient accurately. Instead, it will tumble 
when the output of some filtering of the chemoreceptor histories crosses some threshold 
which indicates that it is headed sufficiently in the wrong direction. The problem of the 
optimal strategy is to determine the filter, F(t) and the threshold. This problem will be 
soluble for two different assumptions about the nature of the tumbles. First, we will treat 
the case of completely disorienting tumbles of finite duration (r tum bic), as for the high signal 
to noise ratio case, and then the case of instantaneous tumbles which do not completely 
disorient the bacteria. 

Let us introduce some notation. In the low signal to noise ratio limit, the filtering should 
be linear and we denote the filtered output by z/(t), defined by: 



where c(t) is instantaneous concentration inferred from the fraction of occupied receptors, 
0(t). Defining the concentration at which half of the chemoreceptors will bind attractant as 



C. Low Signal to Noise Ratio Strategy 




(31) 



C1/2, we have: 
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e« - ^|yC 1/2 (32) 

It is convenient to regard c(t) as the sum of two terms: 

c(t) = u(t) • Vc + t)(*) (33) 

and denote their contributions to u(t) by 

r°° i i - 

s(t) = I dt F(t )v(t) ■ Vc (34) 
n 



o 



(f) = rdt'F(t')f,(t) (35) 
Jo 

The ?7 term represents the fluctuations in the chemoreceptor occupancy associated with the 
stochastic binding of attractant molecules to receptors. 

Now let us compute the "return" , (v ■ Vc) to be maximized. For this, we anticipate the 
result that the characteristic time constant of the filtered "noise", n(t), is short compared 
to that of the filtered "signal", s(t) (see Appendix A). We take the "noise" to be Gaussian 
since it is a linearly filtered version of the noise in the chemoreceptors, and is expected to 
be approximately Gaussian via the Central Limit Theorem (there are many receptors each 
making small contributions to the total noise). 

In this case, the rate at which the filtered output, f(t), crosses the threshold, T, is given 

by: 

r « = & IP{JJ 2W 1) (36) 

with r n defined by 



W (37) 

The trajectory of a bacterium emerging from a tumble at time t = is therefore weighted 
by: 

W (t,{ S }) = e- / ^ t ' r ( s(t ' ) ) (38) 
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The return on the strategy is then given by: 

V ' r turahlc + (Jdt'w(t',{s})) V ; 

where the numerator is the mean distance swum up the gradient of attractant in a run and 
the denominator is the mean time required for a run. 

Let us begin with the computation of the denominator. The typical length of a trajectory 
is given by: 



(jdtw(t,{s})) (40) 



This can be computed to lowest order in the signal to noise ratio more or less straight- 
forwardly; the only caveat is that while s(t) is formally proportional to the gradient of 
concentration, there is the possibility that the coefficient of proportionality could diverge if 
the filter, F, has a long tail. In fact, s(t) can be partitioned into two pieces: one coming 
from the contribution to the integral from times in the not too distant past, i.e. not much 
longer ago than the rotational diffusion time (J ~ r dt' F(t')v(t — t')-Vc, where r = (2D rot )^ 1 ) 
and another piece coming from the remainder of the integral (x = dt' F(t')v(t — t') ■ Vc. 
Only the second piece can avoid being small for small concentration gradient, and only if 
F(t) decays slowly at long times. In this case, this term yields a contribution to s(t) which 
is effectively static and is Gaussianly distributed since it is essentially the sum of many in- 
dependent contributions from s(t) stretching far into the past (recall that v(t) decorrelates 
on a time scale that is of order (2D rot ) _1 ). The mean value of this contribution, (x), is zero 
and 

(z?) = Z™5 ( j dt I dt'F(t)F(t')v(t - t) ■ Vc v(t - t") ■ Vc) (41) 

Trun ^tumble J J 

= ^ / dtF\t ) (42) 

7~ran ^tumble ^ 

where we have assumed: 

(v(t)v(O)) oc e~ 7|t| (43) 
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and low signal to noise ratio. For the case where each tumble completely disorients the 
bacterium, 7 is given at zeroeth order in the concentration gradient by: 

7 = r'i + 2D TOt (44) 

This second contribution, x, to s(t) is can be thought of as a contribution to the noise 
since it is uncorrelated with the present orientation of the bacterium. In fact, since it and 
the occupancy noise are both Gaussian, it can be integrated out. Let us therefore switch to 
using s to refer only to the first part since that is the part which actually does contain the 
signal. Then we find that: 

(r( S ))™„, ltl part o, , = /^iL=|-exp(-^)exp(-<^i^-) (45) 

= f w> (46) 

where n' is a new effective noise whose statistics are those of the original noise, n, plus a 
static contribution, x. Thus 

(n' 2 > = (n 2 ) + (x 2 ) (47) 



( "' 2> (48) 



(n 2 ) 



\ (r 

Hereafter, we switch to using n to refer to the effective noise rather than using n'. We are 
now also in a position to state our definition of the signal to noise ratio: 

SNR=^^ (49) 

The result is implicitly a function of the filtering scheme and the tumbling rate, however, 
one generally expects the signal to noise ratio to be proportional to v 2 |Vc| 2 / (7777). In fact, 
anticipating our results for the optimal filter and tumbling rate we find that in the low signal 
to noise ratio limit: 

SNR = t; 2 |Vc| 2 (l + fr tumblc ) 

N(8D TOt + 4D rot rr tumb ie -r- 3r 2 r tumb i c )(2D rot + r) 2 

where 



N 

(V(t)v(t )) = 5 exp(-\t - t |/r bare ) (51) 

~ N5{t - t) (52) 

We now consider the return at lowest order in the signal to noise ratio. To zeroeth order 
in the concentration gradient, r is constant over a run and given by: 

f = T &^-m ] (53) 

and the mean run length is just f _1 and the denominator of the return expression, Eq. 
is just r tumble + f- Y . 

The numerator of Eq. 39, can also be computed to lowest order in the signal to noise 
ratio. The redefinition of the "signal" and the effective noise that we used for calculating 
the mean run length is appropriate here also so we continue to use s to refer only to the 
contribution from the recent past. We expand the w appearing in the numerator to lowest 
order in s to find: 

(J dtv{t')-Vc w(t',{s})) ~ y2 \ Vc \ 2 dt ' ( cosQ ( t ') q^ c ^ ex P f~ Jl r (t")dt"^j) (54) 

^2| Vc | 2 |r r dt ' e -rt' [* dt " dt"F(t")(cose(t)cos(e(t" -t")) (55) 
oT Jo Jo Jo ' v ' 



-2> 2 |Vc' 2 



where we have used 



3(n 2 )(2D Iot + f) Jo 



dr rT 



/ F(t")exp(-(2D mt +f)t") (56) 
Jo 



dT (n 2 ) 

The return is therefore given by: 



(57) 



v 2 \Vc\ 2 T , . /*oo 
= L-L-— / f) / dt F t exp - 2L> rot + f t 59 

To maximize the return we require that 8r ^ rn = and <5 ^" f r ° = 0. 
For the T equation we require 
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1 1 d(n 2 } ( d\nf(r) dt"T F(t"') exp(-(2Aot + f)t"') \ df 

T (n 2 ) dT { 8r f™ dt"'F(t'") exp(- (2 Aot + ?)?") J dT 1 } 



T 2 ( 2Aot - r 2 r tumble f / °° dt'" t'" F (t'" ) exp(-(2£> rot + f)t" 



-i 



(61) 



(n*) V (! + tumble) (2 A-ot + r) J oo C /rF(r)exp(-(2A-ot + f)r) / 
where we have neglected subleading terms in the signal to noise ratio. 

Now we need to set up the filter equation and solve the two simultaneously. For the filter 
equation we need to know: 



Sf 1 



_f_\ ( T 2 - (n 2 ) 5(n 2 } ^(n 2 ) \ 



KT&k (62) 



8F{t) 2\{n 2 ))\ (n 2 ) 6F(t) n 5F(t)J 

We will also need to know and f^y, for which we need to specify the nature of the 
bare noise, r\. We take (rj) = and 

N 

(V{t)v(t )) = ~ exp(-\t - t |/r bare ) (63) 

~ N5(t - *') (64) 

since the bare noise time constant is much shorter than any other time scale in the problem 
(see Appendix A). In this case, 

§t ~ - 2Np{v > + 3(1+f n!2!^ 1+f ) F(f) (65) 

where these equations are valid for regions where F is slowly varying compared to the bare 
noise time and to leading order in the signal to noise ratio. The filter equation is given by 

( dlnfjr) _ £ dffFjf) exp(-(2Aot + f)t'") \ 6f 
{ df !™dt"'F(t"')exp(-(2D Iot +f)t"') ) 5F(t) 

1 S(n 2 ) exp(-(2Aot + r)t) 



(n 2 )5F(t) ' J oo F(r)exp(-(2Aot+f)r') 
exp (-(2D ro t + f)t) 1 6(n 2 ) ( (n 2 )\ 5f 



f °° F(t"') exp (-(2Aot + f)t"') (n 2 ) 6F(t) \ fT 2 J 6F(t) 
First, consider the behavior of the filter for times large compared to (2D rot + f) -1 . In this 
case the left hand side of Eq. ^8| is negligible. Further, the contribution to ^J^y from the 
jpjf] is small, provided that F is slowly varying in this region. So the filter equation requires, 
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(68) 



2v 2 \Vc\ 2 

FW = 3iV(l + fr tumble )(2Aot + r) F(t) (?0) 



The absolute scale of the filter is arbitrary given the definition of T used in Eq. ^0], so we 
may take: 



F(t) ~ exp — „ 



2v 2 \Vc\ 2 



(71) 



\ 3iV(l + fr tumblc )(2Aot + r) 

The filter thus has an extremely long tail determined by the square root of the signal to 
noise ratio. 

The behavior for intermediate times, times of order (2D TOt + f) _1 , is more complicated. 
Now we must retain the exponential term in Eq. |6^, however, the term in |^ coming from 



5F 

is still negligible and we may also neglect the term in that vanishes with the signal 
to noise ratio. In this case we have: 

( 1 | 1 \ ^ 2 ) _ exp(-(2D rot + f)t) 



\2{n 2 ) 2T 2 j5F(t) f™ df F(t"') exp (-(2£> rot + f)t 

(1 1 \ r°° 
+ Y 2 ) 1 rft "' F(t "' ) exp (-( 2D ^t + r)t") (73) 

This requires that 

F(t) = A + Bt + Cexp(-(2D Iot + f)t) (74) 

Matching the filter onto the result for long times requires that A = 1 and B = 0, while 
-F(O) = requires C = — 1 unless F were to vary very rapidly in the short time region, 
which is clearly not optimal because of the additional contribution to the noise that would 



result. In this case, C = — 1 imposes a self-consistency condition. Eq. 73 requires that 



{(n 2 ) + T 2 ) N(2D mt + f) (75) 

where we have assumed that the very short time 
contribution to Jq 00 dt"'F(t"') exp (— (2D rot + f)t"') is negligible. The contribution to (n 2 ) 
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from the the intermediate and long time parts of the filter is M^±r) so 

that there must 

be a short time contribution to (n 2 ) given by; 

/ 2\ - N\2D rot + ff 

\n ; short times - _ 2iy(2Aot + f) \«>) 

The short time behavior of the optimal, low signal to noise ratio filter depends on the 
characteristics of the bare noise. The short time equation for the optimum filter is: 



/ r°° ,,,,, ,„ s , , _ . II 1 \ 5(n 2 ) 9 1 5(h 2 ) , . 

2(/ o *^)exp(-(2 Dl „, + r-) i )) = 55 + 7^ )M + ^ikf) (77) 



T 2 (n 2 ) J SF{t) n T 2 5F{t) 



where 



S(n 2 ) 2 1 r°° 



Tbare r bare 

<5(n 2 ) 2 2 



^ft)- — rdt'F(t')exp(-\t-t'\/r haxe ) (78) 



T bare T b. 



F(t) (79) 



r 4 

bare 



7— / dtF(t ) 6Xp(-|t - t'|/ T bare) 

Jo 



For small r bare , we may neglect the left hand side of Eq. [77|, and the equation is the solved 
for a filter of the form 

F(t) ~ A (1 - exp(-t/r bare )) (80) 
where A must be chosen so that 



1~n ^ba 



J"2 

\ 1 + 7^\ ( 81 ) 



This requires: 



^2 W 2 D „ t+ ,)^&±lL (82) 
This choice of A must also fulfill the self-consistency equation from intermediate times, Eq. 



73| . A straightforward calculation demonstrates that, to leading order in the signal to noise 
ratio, both self-consistency requirements are satisfied by this choice of A in the limit of small 
Tbare- We thus arrive at the final form of the filter: 
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2T 2 - iV(2D rot + f) ; I \ 3iV(l + fr tumblc ) (2 D rot + 



2v 2 \Vc\ 



-2r> f2Z> If) ^(2Aot + f) 
- exp ((2D rot + f)t) 



Notice that, as expected based on Berg and Purcell's original arguments ||11|| , the time 
scale of the filter is set primarily by the rotational diffusion time of the bacterium. Our 
filter is also similar to that found in [111] , where a different criterion for optimality was used. 



However, the long time behavior found here is somewhat different from that of M] and from 
that expected based on the arguments of JTT]. We find that, in the low signal to noise ratio 
limit, a temporal filter extending for significantly longer than the rotational diffusion time 
is optimal. For weak signal strengths some useful information is gained from averaging in 
measurements made much longer ago than the rotational diffusion time scale, and thus the 
optimal filter includes these times. This is not what one would have naively expected based 



on the arguments of [11]. A similar effect was found in [14], where the tail of of the optimal 
filter was found to extend to infinitely long times, however, that result was obtained in the 
strictly zero signal case, rather than in the limit of small signal as for our filter. We find 
that it holds only for the case of strictly zero signal to noise ratio, (either for our definition 



of optimality or that used in [T|]) and that, for finite signal to noise ratio, the optimal filter 



involves an additional time scale, y3iV(l + fr tum bic)(2-D r ot + f)/(2v 2 \ Vc| 2 ), depending on 
the rotational diffusion time and the signal to noise ratio. The tail of the filter decays on 
this time scale, which diverges with the inverse of the square root of the signal to noise ratio 
but is of the same order as the rotational diffusion time for finite signal strengths (where the 
calculations considered here are not strictly valid but should be qualitatively correct). Our 
results therefore unify the conclusions of of fTl] and [14| : for moderate signal strengths, Berg 



and Purell's argument that the filtering time scale will be of order the rotational diffusion 
time will be correct, but at low signal to noise ratio a new, extremely long, averaging time 
scale is optimal. 
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To compare Eq. |83| to experimental results on the filtering strategies used by the bacteria 
we require the equations which determine f,T, (n 2 ) and r n . These are given by Eq. |53|, 



T = J(n*) 2 ( 1 +^ umble )(2£) rot + f) 
V 4D rot - 3r - 5r 2 r tumb ie 



^ - - r l (85) 



and Eq. 1ST], respectively. We can combine Eq. Rl and 183] to obtain 



- r N , 8 Aot + 4D rot rrtumbie -r- 3r 2 r tumb ie 

T= VT 2jDrot + r In ^ — ^2 86 

V 2 4 Aot - 3r - 5HT tum ble 



2 , N 8D ro t + 4Z} ro trrt um bio — r — 3r 2 r tum bi e , _^ 

n ) = ~r rr— ( 8T ) 



4 1 + rr tumbl , 

The remaining equations must be solved numerically and for the typical values ||: 

numbic ~ -15 sec (88) 

2D mt ~ .3 sec" 1 (89) 

r bare ~ 1msec (90) 

the resulting f is: 

f ~ .17sec _1 (91) 

This value depends only very weakly on the poorly known Tb ar e and the experimentally 
measurable r tum bi e f24|. It depends strongly on D rot and also on the assumption that each 



tumble totally disorients the bacterium. For this reason it is worth considering the problem 
of tumbles which are not perfectly disorienting. 

The low signal to noise ratio is also solvable for tumbles which do not perfectly disorient 
the bacterium in the limit of vanishing T tum bi e - If we define 

Z = 1 - (Wbefore^aftcr) (92) 

then we find that the return (Eq. |10|) becomes: 
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return = ' ' — — / dtF{t) exp (-(2D rot + f)t) (93) 

In this case, all of the arguments for the case of perfectly disorienting tumbles go through 
as before except that: (1) the return is scaled by z, (2) f must be replaced everywhere by zf 
and (3) r tum bi c is to be set equal to zero everywhere. Roughly, the resulting optimal value of 
f will be 1/z larger. Experimentally, it appears that z~l/2l and for this value we find: 

f ~ .37sec _1 (94) 



D. Comparison with Experiments 

Many of the features we would expect in the behavior of bacteria implementing the op- 
timal strategy are directly comparable to the observations of Ref. || and ||25|| , where the 



behavior of free swimming E. Coli in spatial and temporal gradients of various chemoattrac- 
tants, as well as in the absence of such gradients, was studied. 

In Ref. [[5], free swimming E. Coli in the absence of gradients were found to have a distri- 
bution of run times that was approximately exponential with a time constant of about .85 
sec. The distribution of run times could be made almost perfectly exponential by rescaling 
the run times by the mean run times of the individual bacteria. This form for the distribu- 
tion of run times for an individual bacterium is in agreement with the form we find in the low 
signal to noise ratio limit, the relevant limit in this case. The fact that the time constant is 
different for different bacteria is also a natural for the optimal strategy because the bacteria 
differ in their rotational diffusion constants, and, therefore, different bacteria should choose 
different rates for initiating tumbles. It would be very useful to see if differences in the 
rotational diffusion constants of individual bacteria correlate with their different tumbling 
rates, a question not investigated in ||. Assuming that the disorientation due to tumbling 
remains fixed, the tumbling rate from the optimal low signal to noise ratio strategy should 
be roughly proportional to the rotational diffusion constant, if the bacteria are pursuing the 
optimal, adaptive strategy. 
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Whatever the nature of the correlations between rotational diffusion and rates of tumble 
initiation, the mean rate of tumbling observed, ~ 1.2 sec -1 , is anomalously large. It is a 
factor of 3 larger than the rate we would expect even if we take into account the correlation 
between orientations before and after tumbles reported in ||. It is possible that there is 
significant error in either the value of D rot or of z that we have used and it would be of 
great value to have precise experimental determinations of these from tracking experiments, 
since they are both experimentally directly measurable. However, a factor of 3 appears to 
be too much to be the result of inaccuracies in these values, and Ref. f| quotes a mean 
change in orientation from beginning to end of run of only 23 degrees, implying that the 
bacteria really do run for times significantly shorter than the time which disorients them. 
In the low signal noise ratio limit this is not optimal; however, this behavior may be the 
result of the experiments involving bacteria in the absence of any chemoattractant, rather 
than merely in the absence of gradients. In fact, bacteria in a uniform solution of 10 -4 molar 
serine have an exponential run distribution with a time constant that is roughly three times 
longer than was found in the absence of serine [[J. This agrees rather well with the value 
expected for the optimal strategy, however, it should be noted that a uniform concentration 
of aspartate, a different chemoattractant, was not found to have the same effect. Clearly it 
would be desirable to have more tracking work done in uniform or nearly uniform solutions 
of chemoattractant that are as similar as possible to the natural environment of E. Coli in 
order to settle this question. This would enable us to determine whether the tumbling rate 
is really anomalously large compared to the optimal. 

In addition to determining the distribution of run times in the absence of chemoattrac- 
tant, the observations of || also measured some of the effects of small gradients on this 
distribution. One of the interesting things found was an indication of a peculiar asymmetry 
in the response of the bacteria to small gradients. E. Coli tumbled less often when swimming 
up the gradient of chemoattractant, but not more often when swimming down the gradient 
of chemoattractant. We have seen that at high signal to noise ratio, bacteria following the 
optimal strategy will be surprisingly reluctant to tumble because of the finite amount of 
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time required to tumble and the fact that there is some chance that rotational diffusion will 
improve their prospects. This may explain some of the observed reluctance to tumble if 
the bacteria are in a medium signal to noise ratio regime. In fact, the asymmetry is only 
clearly observed for runs that are longer than 1.5 sec, so perhaps a minimum integration 
time, and the accompanying boost in signal to noise ratio, is required for the asymmetry. 
Any detailed attempt to explain the asymmetry would require more detailed measurements 
of the responses of E. Coli than currently available for free swimming bacteria, and this is 
another area where the collection of more data on free swimming bacteria would be of great 
value. 

At low signal to noise ratio, there is a possibility for asymmetry in responses other than 
the cost of tumbles: we expect an exponential, not linear, dependence of the tumbling rate 
on the sensory input. In fact, Ref. |25| observed the response of free swimming E. Coli to 



temporal gradients of glutamate and found that their results were best fit by an exponential 
dependence on the rate of change of receptor occupancies, which in the concentration region 
treated implied an exponential sensitivity to gradients in concentration. They did not, 
however, claim to have ruled out a linear dependence and further measurements of the 
response of the tumbling rate in bacteria adapted to a low signal to noise ratio environment 
would be of great value. It would also be of great value if bacteria could be placed in 
a range of spatial gradients of chemoattractant and then stimulated with an additional 
temporal gradient to see if the response crossed over from one appropriate for the low signal 
to noise ratio case to one appropriate for the high signal to noise ratio case, i.e. tumbling 
at fixed value of the estimated angle from the concentration gradient. 

The predictions we have made for the optimal filter to be used in the low signal to 
noise ratio limit can be compared to the results of [|26] , [T0| , [T2]1 , where tethered |20| bacteria 



were exposed to impulse like bursts of chemoattractant. If the strategy the bacteria employ 
involves linearly filtering some function of sensory inputs, then, the response a time t later 
under these conditions is related to E(t), and so the temporal properties of the derivative 
of the filter can be determined from these experiments. Evidence for the very long tail in 
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the derivative of the filter expected of the optimal strategy for low signal to noise was not 
found in these experiments. However, the filter found there does extend about 4 sec into 
the past, which is significantly longer than the time scale for the bacterium to disorient. 
The improvement in chemotactic performance that would result from a longer integration 
time is very small, while the difficulty of building a faithful, long term memory is clearly 
substantial, so the result is not surprising. The long time behavior of the filter may still be 
adaptive but the limit of truly small signal to noise ratio appears to be beyond the range (if 
any) of that adaptation. 

At intermediate times the derivative of the optimal filter contains a term with an ex- 
ponential decay rate given by 2D TOt + zf <~ 1 sec -1 (for the realistic assumptions of short 
tumble times and finite disorientation during tumbles). A feature of almost exactly those 
characteristics is seen in the filter found in ||10|| , and it appears that the optimal filtering 
strategy may in fact be rather close to the filtering strategy inferred from the response ob- 
served in those experiments. However, one should be cautious since the agreement is the 
result of using the experimental value for f which is not in agreement with theoretical expec- 
tations (unless our value of D mt or z is badly off), and there are worries about the adaptive 
state of tethered bacteria []27|] . Also, the tethered experiments were interpreted in terms 
of a linear, but thresholded, response in the rate of transitions from running to tumbling 
and vice versa; this is not readily reconciled with the proposed deterministic strategy and 
the resulting exponential dependence of the effective rate of such transitions on the filtered 
signal. 

It is probably not sensible to compare the behavior of the filter at short times (100 msec or 
shorter) with theoretical expectations since responses on this time scale are not particularly 
important for chemotactic performance, while constraints due to the actual physical signal 
processing mechanisms of the bacteria are severe on this time scale. 

In summary, the filter observed has behavior very similar to what is expected for the 
optimal filter at time scales of order a second. On longer time scales it decays, faster than 
optimal rate for low signal to noise ratio, but at a rate that is clearly slower than the natural 
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time scale for the bacterium to disorient, in qualitative agreement with that characteristic of 
the optimal filter. There is some question as to the adaptation state of the tethered bacteria, 
but since they spend most of their time in the complete absence of chemoattractant signals, 
the most reasonable proposal is that they are adapted to low signal to noise ratio. The 
tethered experiments do not make any attempt to identify adaptation to signal to noise 
ratio or diffusion constant in the characteristics of the filter. In fact, they are ill suited 
to such a test since the tethered bacterium experiences such unnatural conditions that its 
adaptive state is difficult to determine. 

On the other hand, some experiments suitable for testing the form of the filter on free 
swimming bacteria have recently been performed using photoreleased chemoattractant and 
repellant . In those experiments, bacteria swimming freely in the absence of spatial gra- 



dients were exposed to step like changes in concentration by photoreleasing caged chemoat- 
tractants and repellants at the location of the bacteria. This form of stimulation provides 
a delta function in the derivative of the concentration so that the response a time t later 
is a measure of the filter, F(t), used in the low signal to noise ratio environment. These 
experiments find a memory time for the system slightly longer than that of the tethered 
experiments, around 5 sec. This is marginally closer to the long time tail is expected for 
the optimal filter. On the other hand, there is no indication of the intermediate time (~ 1 
sec) feature found in the data of [26|, although the data in the case of [2£| are somewhat 



noisy and the applied stimulus is so large, inducing a tumble in 100 msec with nearly unit 
probability, that a feature with the expected rise time of ~ 1 sec is probably not excluded. 

It would be very useful to have more photorelease-based studies of the response of free 
swimming bacteria aimed specifically at measuring the properties of the filter on time scales 
of a fraction of a second to several seconds. One would like to know the filter properties as 
a function of the spatial gradients to which the bacteria are pre-adapted, and, if possible, 
one would also like to look for correlations between the time constants of the filter used, 
the rotational diffusion constant of the free swimming bacteria and their distribution of 
run times. Experiments of this type could answer the important open questions of the 
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adaptability of the form of the filter and the tumbling criterion directly. This would be a very 
important step in determining the extent to which E. Coli achieves the optimal chemotactic 
strategy. We believe that this would provide interesting and important information about 
the limits of the sophistication of sensory processing in single celled organisms. 
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I. APPENDIX A 



The correlation time for the chemoreceptors is expected to be of order 10~ 3 sec or shorter. 
This time scale is results from several considerations. First, the typical binding time for a 
chemoreceptor is about 10 -4 sec [11]. This will set the correlation time for the chemoreceptor 
inputs if: (1) it is a much larger time than the time for attractant to diffuse away from a 
receptor, and, (2) there is relatively little correlation between receptors. To see that (1) 
is satisfied we approximate the time to diffuse away from the receptor (that is the time to 
diffuse far enough away to have small chance of recapture) as the Slze of the receptor 2 w hich 

° J ' diffusion constant of attractant 

is typically 10~ 14 cm 2 /10 _5 cm sec -1 ~ 10~ 9 sec, vastly smaller than the rotational diffusion 
time, even if we have underestimated the distance out to which recapture is important by a 
sizable factor. Individual receptors therefore decorrelate on a time scale given by the binding 
time. What about the ensemble of receptors? 

To see that the ensemble of receptors decorrelates on a time scale comparable to the 



binding time, recall the arguments of |Tl] regarding the rate of capture of diffusing molecules 
by a large number of small perfectly absorbing sites on the surface of an impermeable sphere. 
The inbound current for N patches of linear size s on a sphere of size a is given by: 

Ns 

J = AnD Coo a-— (95) 

Ns + na 



where is the concentration of signaling chemical per cubic centimeter at infinity. In |TT 



it was emphasized that this differs from the current to a perfectly absorbing sphere of radius 
a only by a factor of N ^ na , which can approach one for reasonable choices of N, s and a. 
Here, we note that this occurs because the current ^TiDc^a AT N , S differs from that for N 
independent disks only by a factor of N ^ 7Ta , which is about one-half for the parameters of 
||11|| . Clearly, the reason for the reduction from the value for N independent receptors is 
that some of the molecules absorbed by a given receptor would have contacted others in its 
absence. In fact, the average molecule that contacts an absorbing site would make contact 
with another binding site with probability N ^ na ■ For receptors which bind and release 
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attractant, this implies that the average molecule which binds to one receptor will bind to 
V(l ~~ Ari+ S 7ra ) = NS wa a ~ ^ others, so that, defining the number of receptors with attractor 
bound at time t to be x(t) and assuming Poisson statistics: 

[ T dt(x(t)x(0)) - T(x) 2 = ( Ns + na \ Th . md{x) (96 ) 
Jo \ na / 

~ 2 r bind (x) (97) 

If the receptors decorrelated on a time scale, r decorr , which was much longer than r bind , then 

f T dt{x(t)x(0)) - T(x) 2 ~ e Trft exp(-t/r dccorr ) (x(r bind )a;(0)> (98) 

~ er dccorr (x(r bind )x(0)) (99) 

implying 

(x(T bind )x(0)> ~ (100) 

C^decorr 

< 1 (101) 

and the correlations among the receptors must be weak if their correlation time is much 
larger than the binding time. Conversely, if r dcC orr ~ Tbind, then the correlations need not 
be weak, but the two time scales are comparable so that r bind is still the appropriate time 
scale. 

In practice, the time to diffuse far enough away from all receptors should be roughly 
the size of the bacterium 2 which ■ typically io- 8 C m 2 /10" 5 cm sec" 1 ~ 10- 3 sec or ten times as 

dmusion constant ot attractant j r j i 

long as the receptor binding time. In this case, most of the correlations among chemoreceptor 
inputs have decayed by .1 msec but some weak correlation persist out to 1 msec. Both times 
are still very short compared to the relevant rotational diffusion time scale, and we can 
approximate the noise in the binding of the chemoreceptors to be white on the other time 
scales of interest. 
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FIGURES 




tumble 

FIG. 1. Plot of the optimal angle at which the bacterium should initiate tumbles versus the 
product of its rotational diffusion constant and the duration of the disorienting tumbles. 
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